Observation of a robust peak in the glue function of the high-Tc cuprates in the 50-60 

meV range. 



E. van Heumen/ E. Muhlethaler/ A.B. Kuzmenko/ D. van der 
Marel,^ H. Eisaki,^ M. Greven,^ W. Meevasana,^ and Z.X. Shen^ 

' Departement de Physique de la Matiere Condensee, Universite de Geneve, 
quai Ernest- Ansermet 24, CHI 211 , Geneve 4, Switzerland 
Nanoelectronics Research Institute, National Institute of Advanced Industrial Science and Technology, Tsukuba, Japan 
'^Department of Physics, Applied Physics, and Stanford Synchrotron 
Radiation Laboratory, Stanford University, Stanford, OA 94305 

We take advantage of the connection between the free carrier optical conductivity and the glue 
function in the normal state, to reconstruct from the infrared optical conductivity the glue-spectrum 
of ten different high-Tc cuprates revealing a robust peak in the 50-60 meV range and a broad 
continuum at higher energies for all measured charge carrier concentrations and temperatures up 
to 290 K. We observe an intriguing correlation between the doping trend of the experimental glue 
spectra and the critical temperature. 



The theoretical approaches to the high Tc pairing 
mechanism in the cuprates are divided in two main 
groups: According to the first school electrons form pairs 
due to a retarded attractive interaction mediated by vir- 
tual bosonic excitations in the solid |TJ |2 [31 [S]. These 
bosons can be lattice vibrations, fluctuations of spin- 
polarization, electric polarization or charge density. The 
second school concentrates on a pairing-mechanism en- 
tirely due to the non-retarded Coulomb interaction [5] or 
so-called Mottness[7 . Indeed, optical experiments have 
found indications for mixing of high and low energy de- 
grees of freedom when the sample enters into the super- 
conducting state H Uni HB . 

An indication that both mechanisms are present was 
obtained by Maier, Poilblanc and Scalapinofl^, who 
showed that the 'anomalous' selfenergy associated with 
the pairing has a small but finite contribution extend- 
ing to an energy as high as U, demonstrating that the 
pairing-interaction is, in part, non-retarded. They ar- 
rived at the conclusion, that "for the cuprate materials, 
the relative weight of the retarded and nonretarded inter- 
action remains an open question. Thus, the continuing 
experimental search for a pairing glue in the cuprates is 
important and will play an essential role in determining 
the origin of the high-Tc pairing interaction." Aforemen- 
tioned glue can be expressed with the help of the density 
of states of these bosons multiplied by the electron-boson 
coupling, a'^F(uj) for phonons and I'^xi^) for spinfluctu- 
ations, in this Letter represented by the general symbol 
n([j) . An important consequence of the electron-boson 
coupling is, that the energy of the quasi-particles relative 
to the Fermi level, i^, is renormalized, and their lifetime 
becomes limited by inelastic decay processes involving 
the emission of bosons. The corresponding energy shift 
and the inverse lifetime, i.e. the real and imaginary parts 
of the self-energy, are expressed as the convolution of the 
'glue-function' n(w) with a kernel K{£^,uj,T) describing 



the thermal excitations of the glue and the electrons [l3] 



S(0 = J K{^,LU,T)JI{Lj)dLJ 



(1) 



In the absence of a glue and of scattering off impu- 
rities the effect of applying an AC electric field to the 
electron gas is to induce a purely reactive current re- 
sponse, characterized by the imaginary optical conduc- 
tivity 4TTa{uj) = iujp/u!, where the plasma frequency, ojp, 
is given by the (partial) f-sum rule for the conduction 
electrons. The effect of coupling the electrons to bosonic 
excitations is revealed by a finite, frequency dependent 
dissipation, which can be understood as arising from pro- 
cesses whereby a photon is absorbed by the simultane- 
ous creation of an electron-hole pair and a boson. As 
a result, the expression for the optical conductivity in 
the normal state, 47rcr(tj) = iUp/luj + 2Y.opt{i^)}, now 
contains a memory function|14j equivalent to an 'optical 
self-energy'. A particularly useful aspect of this repre- 
sentation is that Sopt(a') follows in a straightforward way 
from the experimental optical conductivity. The optical 
self-energy is related to the single particle self-energies 
by the expression [T5] 



2Sopt(a;) 



+ s*(0-s(C + c^)' 



rfQ -1 (2) 



The central assumption in the above is the validity of 
the Landau Fermi-liquid picture for the normal state. 
The aforementioned strong coupling analysis is there- 
fore expected to work best on the overdoped side of the 
cuprate phase diagram, where the state of matter appears 
to become increasingly Fermi liquid like. If antiferromag- 
netism is necessary to obtain the insulating state in the 
undoped parent compounds, as has been argued based on 
the doping trends of the Drude spectral weight [TB], the 
strong coupling analysis may in principle be relevant for 
the entire doping range studied. However, in the limit 
of strong interactions aforementioned formalism needs to 
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be extended, e.g. with vertex corrections, and eventually 
breaks down. We therefore define the function n(a;) as 
the effective quantity which, in combination with Eqs. 
[l]and[2] returns the exact value of Eopt(a;) for each fre- 
quency. Defined in this way n(w) captures all correlation 
effects regardless whether the system is a Fermi-liquid or 
not. This becomes increasingly relevant when the doping 
is lowered below optimal doping. 

Here we take advantage of the connection between the 
temperature and frequency dependent conductivity in 
the normal state and the glue-spectrum to test experi- 
mentally the consequences of the standard approach, to 
check the internal consistency of it, and to determine 
the range of doping where internal consistency is ob- 
tained. For a d-wave superconductor, the momentum 
dependence is essential to understand the details of the 
pairing. This, of course, is difficult to handle for opti- 
cal spectroscopy which is inherently a momentum inte- 
grated probe. Nevertheless, optical spectra provide the 
important information on the energy scale of the bosons 
involved and on the doping and temperature evolution. 
We use a standard least squares routine to fit a histogram 
representation of n(a;) to our experimental infrared spec- 
tra. The quantity n(ix') is shown in Fig. [ijfor optimally 
doped HgBa2Cu04+5 (Hg-1201)[17 for T = 290 K, to- 
gether with the optical self energies calculated from this 
function at three different temperatures. For 290 K the 
theoretical curve runs through the data points, refiecting 
the full convergence of the numerical fitting routine. It 
is interesting to notice, that the shoulder in the exper- 
imental data at 100 K is reproduced by the same 
function as the one used to fit the 290 K data. In other 
words, the strong temperature dependence of the exper- 
imental optical spectra is entirely due to the Fermi and 
Bose factors of Eqs. [l] and |2] It can be excluded that 
the shoulder at 80 meV is due to the pseudo-gap, since 
a gap is certainly absent for temperatures as high as 290 
K. The shoulder is therefore entirely due to coupling of 
the electrons to a mode at approximately 60 meV. On 
the other hand, the considerable sharpening of this fea- 
ture for temperatures lower than 100 K finds a natural 
explanation in the opening of a gap, as illustrated in the 
inset of Fig. [T] This example confirms the close corre- 
spondence between the features in Sopt(w) and in n(a;) 
pointed out in Ref. |18j . In particular the broad maxi- 
mum in Sopt (uj) has its counterpart in the high intensity 
region of n([x)) terminating at 290 meV. 

As summarized in Fig. [2] we have analyzed previously 
published optical spectra of 6 different samples belonging 
to different families of materials, i.e. optimally doped 
Hg-1201[Tr and BizSrsCaaCuaOio+a (Bi-2223)[II], as 
well as four BisSraCaCuaOg+a (Bi-2212) crystals [9l [10] 
with different hole concentrations. In addition, we ana- 
lyzed new data for four Bi2Sr2Cu206+5 (Bi-2201) crys- 
tals with different hole concentrations . 
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FIG. 1: Experimental optical self energy of HgBa2Cu04+5 
for 3 selected temperatures (open circles). The solid curve 
at 290 Kelvin is obtained from a fit of n(tj) , shown as the 
dashed surface. The sofid curves at 100 K and 20 K were 
calcufated with the same n(w) function corresponding to 290 
Kefvin. This proves that the seif energy feature between 80 
and 100 meV (a shoufder at 100 K and a peak at 20 K) is 
caused by the prominent peak in tl{u) at approximatefy 60 
meV. The sharpening of this feature at iow temperature is due 
to the superconducting gap, an aspect not captured by Eq. [2] 
and therefore not reproduced in the cafcufated solid curves. 
In the inset the gap-induced sharpening is iifustrated by the 
opticai seff energy without (bfack) and with (red) a 15 meV 
superconducting gap, calculated using Ailen's retation[T5]. 



Excellent fits were obtained for all temperatures, but 
the n(a-') spectra exhibit a significant temperature depen- 
dence, in particular at the low frequency side of the n(a;) 
spectrum. Since all thermal factors contained in Eqs. 
[T] and [2] are, in principle, folded out by our procedure, 
the remaining temperature dependence of Il(a;) reflects 
the thermal properties of the 'glue-function' itself. Such 
temperature dependence is a direct consequence of the 
peculiar DC and far infrared conductivity, in particular 
the T-linear DC resistivity and lu/T scaling of Ta{uj,T) 
at optimal doping[20^. For the highest doping levels both 
n(Lj) and its temperature dependence diminish, which is 
an indication that a Fermi liquid regime is approached. 
The most strongly underdoped sample, Bi-2201-UDO, ex- 
hibits an upturn of the imaginary part of the experimen- 
tal optical self-energy for u; ^ 0. This aspect of the data 
can not be reproduced by the strong coupling expression, 
resulting in an artificial and unphysical peak at w of 
the fitted n(w) function. 

We observe two main features in the glue-function: A 
robust peak at 50-60 meV and a broad continuum. The 
upper limit of n(w) is situated around approximately 300 
meV for optimally doped single layer Hgl201, and for 
the bilayer and trilayer samples. Interestingly the con- 
tinuum extends to the highest energies (550 meV for the 
single-layer samples and 400 meV for the bilayer) for the 
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FIG. 2: Electron-boson coupling function ti-{Lo) for Bi-2201 at four different charge carrier concentrations (10 K, 100 K, 290 K), 
Bi-2212 at four charge carrier concentrations, and optimally doped Bi-2223 and Hg-1201 (100 K, 200 K, 290 K). The dotted 
curve in the lower right panel represents the spin-fluctuation model[51 118]. 
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TABLE I: Strong coupling parameters of the ten compounds. 
The hole-doping is indicated on the first line. From left to 
right: Bi2201 (columns 1 to 4), Bi2212 (columns 5 to 8), 
Bi2223 (8th column) and Hgl201 (column 10). The definition 
of Lj is ln(£j;) = 2X~^ uj~^'!1{lli) ln{uj)dijj. 



weakly overdoped samples, whereas the continuum of the 
strongly doped bilayer sample extends to only 300 meV. 
There is also a clear trend of a contraction of the con- 
tinuum to lower energies when the carrier concentration 
is reduced. Taken together this doping trend suggests a 
possible relation of the glue-function to a quantum crit- 
ical point at a doping concentration slightly higher than 
optimal doping [2(T, ^TP. Hence, part of the glue function 
has an energy well above the upper limit of the phonon 



frequencies in the cuprates 100 meV). Consequently 
the high energy part of n(a;) reflects in one way or an- 
other the strong coupling between the electrons them- 
selves. The coupling constant is obtained from the re- 
lation A = 2 J^Il{uj)/ujdLL>. The result (see Table |l]) 
shows a strong and systematic increase of A for decreasing 
hole concentration, which probably requires a theoretical 
treatment beyond the strong coupling expansion outlined 
in the introduction. 

The most prominent feature, present in all spectra re- 
produced in Fig. [2| is a peak corresponding to an average 
frequency of 55 meV. Perhaps the most striking aspect 
of this peak is the fact that its energy is practically inde- 
pendent of temperature (up to room temperature) and 
sample composition. Moreover, the intensity and width 
are essentially temperature independent. While our re- 
sults confirm by and large the observations of Hwang 
et al. in the pseudo-gap phase [22l [23], the persistence 
of the 50-60 meV peak to room temperature has not 
been reported before. This peak arises most likely from 
the same boson that is responsible for the 'kink' seen 
in angle resolved photoemission (ARPES) experiments 
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FIG. 3: Experimental critical temperature and Tc's calcu- 
lated from the McMillan formula using the experimentally 
measured ri(<^) of Fig. |2|at 290 Kelvin as input parameters. 



along the nodal direction in k-space at approximately 
the same energy [551 US]- The peak-dip-hump struc- 
ture in the tunneling spectra (STS)[171 HHJ HSj has also 
been reported at approximately the same energy. 

It is interesting to correlate our experimentally ob- 
tained glue-function with Tc. Millis, Varma and 
SachdevfSn] have shown, that, if n(cj) can be decomposed 
in an s-wave and a d-wave channel and /i* is considered 
negligibly small for the d-wave channel, then the critical 
temperature of the d-wave superconducting phase tran- 
sition is Tc < 0.83a;exp(— (1 -I- A)/A) where ui is the loga- 
rithmic average of the frequency in n(u;) . We calculated 
this upper limit using the experimental values indicated 
in Table m The T^'s with this formula are in the 100-200 
K range for all samples studied here, and, as shown in 
Fig. |3] they correlate with the experimentally observed 
doping trends of Tc for the single-layer and bilayer high- 
Tc materials. The calculated critical temperatures are 
about twice as large as the experimentally observed ones 
for the Bi2212 series, whereas in the Bi2201 single layer 
compounds Tc appears to be strongly suppressed as com- 
pared to the calculated value. The low value of Tc in 
single layer Bi2201 correlates with a relatively high in- 
tensity of Il{uj) below the 50-60 meV peak, and a very 
low intensity or absence of the 0.25 eV peak. In contrast, 
the three samples with the highest critical temperature, 
Bi2212 OpD88, Bi-2223 and Hg-1201, show the most pro- 
nounced weight at ~ 0.25 eV. 

In summary, the H{io) spectrum obtained from the op- 
tical spectra of 10 different compounds using a strong 
coupling analysis, is observed to consist of two features: 
(i) a robust peak in the range of 50 to 60 meV and (ii) 
a doping dependent continuum extending to 0.3 eV for 
the samples with the highest T^. Finally, we observe an 
intriguing correlation between the doping trend of the 
experimental glue spectra and the critical temperature. 
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APPENDIX 

The inversion of Eq's 1 and 2 allows to extract Il{uj) 
from experimental data of the optical conductivity, or 
related optical spectra. The accuracy of the resulting 
n(a;) spectrum is in practice limited by the convolution 
with thermal factors expressed by Eq's 1 and 2 [3T]. Mi- 
croscopic models giving roughly the same n(w) spectra, 
which differ however in the details of the frequency de- 
pendence of this quantity, may therefor provide fits to the 
directly measured optical quantities, such as infrared re- 
flectance spectra, which at first glance look satisfactory, 
but the remaining discrepancies with the experimental 
spectra may nevertheless be of significant importance for 
the proper understanding of the optical data. It is there- 
fore of crucial importance to test the 'robustness' of each 
fit with regard to the spectral shape of the n(w) function 
imposed by such models. This robustness can be tested 
by including in the fit-routine one or several 'oscillators' 
superimposed on the model function. When the model 
glue function provides a complete description of the elec- 
tronic structure, adding extra oscillators will not result in 
an improvement of the quality of the fit. We have used 
this approach to test functional forms commonly used 
in the literature, in particular the marginal Fermi liquid 
(MFL) model [2j and the MiUis-Monien-Pines (MMP) 
representation of the spin fluctuation spectrumlH]. We 
found that neither of these functional forms describe 
completely the experimental data. 

In search of a more flexible form of II(cj) we used a 
superposition of lorentzian oscillators and found that it 
could be used to describe all available experimental data 
in a consistent manner. The resulting H{io) functions and 
trends are equivalent to those in Fig. [2j From these ini- 
tial tests we concluded that due to the thermal smearing 
expressed by Eq's 1 and 2 our II(a;) spectra can only be 
determined with limited resolution. This lead us to the 
use of a histogram representation, where each block in 
the histogram represents a likelihood to find coupling to 
a mode with a well determined coupling strength. Only 
for the lowest frequency interval (0 < w < wi) a triangu- 
lar shape was used instead of a block, which is necessary 
to avoid problems with the convergence of the integral 
A = 2 Jg uj^^Il{uj)du!. In practice the output generated 
by the fitting routine has low intensity in this first inter- 
val, and the triangles are therefore difficult to distinguish 
in Fig. |2] 
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To give an example: The block centered at 55 meV 
seen in the Hg-1201 sample in Fig. [2] has A ~ 1 and 
a width of about 30 meV. Our histogram representation 
implies the presence of a coupling to one or several modes 
between 45 meV and 75 meV with an integrated coupling 
strength of 1. The histograms thus constitute the most 
detailed representation of n(a;) given the precision of our 
experimental reflectivity and ellipsometry spectra. 

Examples of experimental reflectivity data together 
with the fits are shown in Fig. |4]for a selection of repre- 
sentative data sets spanning the entire doping and tem- 
perature range. As the fitted curves are within the limits 
of the experimental noise, further reduction of x^i while 
in principle possible by fitting the statistical noise of the 
data, can not improve the accuracy of the n(ti;) functions. 
Starting from a Ii{i-u) function we can calculate the opti- 
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FIG. 4: Experimental reflectivity (red and blue lines) and flt 
curves (black lines) for selected samples and temperatures, 
(a): underdoped non-superconducting Bi-2201 (UDO). (b): 
optimally doped Hg-1201 with Tc ^ 97 K (OpD97) IIZI (c): 
overdoped non-superconducting Bi-2201 (ODO). Weak sharp 
peaks, particularly visible for the strongly underdoped sample 
in panel (a) are due to transverse optical phonons, which we 
do not intend to flt. 

cal conductivity, which in turn is fed into standard Fres- 



nel expressions to calculate the experimentally measured 
quantities, i.e. refiectivity and ellipsometric parameters. 
The fitting routine is based on the Levenberg-Marquardt 
algorithm and uses analytical expressions for the partial 
derivatives of the reflectivity coefficient R, and the ellip- 
sometric parameters ip and A relative to the parameters 
describing the n(w) function. The algorithm is based 
around minimizing a functional which is given by. 



(3) 



where R^cui) is an experimentally measured datapoint, 
f{uji,pi, ...,Pn) is the calculated value in this point based 
on parameters pi, ...,pn and the difference between these 
two is weighed by the errorbar ai determined for Riuji). 
For a given set of reflectivity and ellipsometry data at one 
particular temperature, using a standard PC, the itera- 
tion takes about 3 hours until convergence is reached. Al- 
though the Levenberg-Marquardt least squares method is 
an extremely powerful method to find the minimum of 
in a multidimensional parameter space it also has diffi- 
culties with finding the best solution to the non-linear 
problem Eq's 1 and 2. For each individual sample and 
temperature displayed in Fig. [5] several tests have been 
performed where in each test the optimization process 
was started from a different set of starting parameters to 
ensure that has converged to the global minimum in 
parameter space. To give some idea of the robustness of 
our method we will here discuss one representative ex- 
ample: optimally doped Hg-1201. 

The models are evaluated based on the minimum found 
for . A comparison of Fig. [s] (a-d) shows that the 
MMP model describes better the optical data then the 
MFL model but that they give similar results if we add 
an extra oscillator to these models. Panels ^5^-f show 
the model independent results mentioned above and are 
very similar to the modified MMP and MFL model. The 
models in these last two panels have the same and 
the comparison in Fig. |5^ shows that the histogram rep- 
resentation realistically expresses the uncertainty in the 
position of the low energy peak, while the correspondence 
between the features in both models remains excellent. 
It is interesting that the model with two oscillators is de- 
scribed by 6 parameters, while the histogram represen- 
tation uses 12 parameters. The fact that the fit-routine 
adjusts the latter 12 parameters in such a manner as to 
produce in essence the two oscillator lineshape, proves 
that the features represented in the righthand panel of 
Fig. [5] are realistic. An extensive discussion of this anal- 
ysis is provided in Ref. [5^ . 
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